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Abstract 

We give a detailed account of the derivation of a master equation for two coupled 
cavities in the presence of dissipation. The analytical solution is presented and phys- 
ical limits of interest are discussed. Firstly we show that the decay rate of initial 
coherent states can be significantly modified if the two cavities have different decay 
rates and are weakly coupled through a wire. Moreover, we show that also deco- 
herence rates can be substantially altered by manipulation of physical parameters. 
Conditions for experimental realizations are discussed. 
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1 Introduction 



The past decade can celebrate having overcome several barriers, the most im- 
portant of them being the objections commonly levelled against the idea of 
the observability of large scale quantum phenomena. The amazing techno- 
logical developments in several areas of physics allowed for the experimental 
verification of essentially quantum effects like tunneling, various kinds of in- 
terference and entanglement phenomena in mesoscopic physics. So we have 
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reached the point where we need to know rather urgently what are the mech- 
anisms which cause decoherence, particularly at the mesoscopic scale, since 
controlling such mechanism has immediate and profound technological con- 
sequences, besides its value for the fundamentals of Quantum Mechanics and 
the (unsolved) problem of the classical limit. 

Modern cavity quantum electrodynamics (cavity QED) sheds light onto the 
most fundamental aspects of coherence and decoherence in quantum mechan- 
ics. One of the main advantages of such experiments is that they can be 
described by elementary models, being at the same time rich enough to re- 
veal intriguing subtleties of the interplay of coherent dynamics with external 
coupling. It represents therefore a unique paradigm for matching theory with 
experiment in the study of quantum coherence. In particular a very impor- 
tant experiment has been proposed [1] and realized [2] where one monitors the 
coherence loss of a superposition state. 

The predominant component of decoherence in cavity QED systems corre- 
sponds simply to the scape or emission of photons, either by absorption into 
the cavity walls or by scattering or transmission into electromagnetic modes 
outside the cavity. This decoherence mechanism fits well into the scheme of 
coupling to extended modes. The description of decoherence through a linear 
coupling to a set of harmonic oscillators seems to work well, both qualitatively 
and quantitatively. 

The purpose of the present work is to investigate dissipation and decoherence 
effects on two field modes by generalizing the usual dissipation mechanism for 
the present situation. The results can be experimentally implemented in two 
distinct ways: in a regime with two modes of one cavity close to resonance with 
an atomic transition, a task already implemented [3], and in an experiment 
involving two cavities, as the one proposed in ref. [4], or in the teleportation 
scheme [5]. 

Our results show interesting physical effects for the situation when the cou- 
pling of each mode with the reservoir is not independent. Despite of being of 
difficult experimental implementation, we show which are the key parameters 
or combinations of parameters which can lead to a decoherence free subspace, 
or, more modestly, to decoherence control. 

Moreover we investigate the coupling of two cavities through a wire and show 
how dissipation can be controlled if the cavities are initially fed with coherent 
states. 

The counterintuitive effects here presented are consequences of quantum de- 
structive interference of coupling constants, a well studied theoretical prevision 
[6] with interesting consequences as electromagnetically induced transparency 
(EIT) [7], lasing without inversion [8] and narrowing of spectral lines [9,10]. 
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The paper is organized as follows. In section II we derive the master equation 
from a Hamiltonian for two coupled oscillators in the presence of an environ- 
ment, under some dynamical hypothesis. Section III presents the exact formal 
solution of the Liouvillian by means of Lie-algebraic techniques. There we dis- 
cuss conditions for dissipation control. In section IV we analyse the case with 
initially coherent states in the cavities and present an interesting effect in the 
case of weak coupling between them. Section V discusses the conditions for 
decoherence control through the study of an initially pure entangled state. 
In section VI we comment on the (difficult) experimental realization of the 
conditions found in the previous sections. 



2 Master equation derivation 



We consider two modes linearly coupled both to each other and to the envi- 
ronment. In both cases applying the Rotating Wave Approximation (RWA), 
what is not too restrictive in the microwave domain [11]. This corresponds to 
the Hamiltonian 

H s = huJaS^SL + hwbb^b + hg (a^b + b*a) , 

k 

Hint = Yl h («fe c l a + a^Cfc) + H h (Afccjb + /3feb f c fc ) , (1) 

k k 

where the operators a (b) and (bt) are annihilation and creation operators 
of the field mode A (B), directly coupled by the last term in eq. (1), where g 
is a real positive coupling constant. The environment is modelled, as usual, by 
a set of harmonic oscillators with creation and annihilation operators c k and 
Cfc, with complex coupling constants to the mode a {b). Much of the 

results to be presented will depend on the existence of correlations between 
«fc and (3k for fixed k, as must become clear later on. 

It is convenient to work with the normal modes of H^: 
H 5 = ft^ia|ai + huj 2 SLlsL 2 , 

2 

H int = Y,Y, h ( a ikCk*i + a*fe a l c fe) , ( 2 ) 
i=l k 

where 
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(u a + uj b ) + J (u a - u b ) 2 + 4g 2 (u a + u b ) - J (u a - u b ) 2 + 4g 2 

"1 = g ' W2 = 2 : 

aifc = ttfc cos + Pk sin tt2fc = — sin + f3 k cos 6*, 
a x = a cos 6 + b sin 0, a 2 = —a sin # + b cos 9, 



cos : 



1 + 



UJ a - u b 



sin 



1 - 



L0 a - L0 b 



y\ yfc a - u b f + v 



( ) 



The operators a^ and a| obey the usual commutation relations for bosons. 

Time evolution of the complete density operator p is given by von Neumann 
equation 



j t p{t) = -^]R int (t),p(t) 



(4) 



where the tilde over the operators tells us that they are written in the inter- 
action picture, 



O ( x i ) { ^ n o t I O exp I H Q t 



(5) 



Integrating eq. (4) and iterating twice leads to 



t f 



p(t)-p (0) = ~f dt'[H mt (if) , p (0)] - ly dt' jdt"[U mt (If) , [H mt (t") , p (0)]] + O 3 (H int ) .(6) 



o o 



Since the high quality factor of the cavity (ies), we consider weak coupling of 
modes A and B to the environment, which allow us to disregard the terms of 
third order in H int , denoted O 3 (ll int ^, in the following calculation. 

Taking the trace over the environment degrees of freedom gives the correspond- 
ing equation for the reduced density operator for the subsystem composed by 
the two relevant modes (this subsystem will be called system) and eq. (6) 
then implies 



1 2 

~ps (t) - ps (o) = -i I dt'Tr E Y: [rj (t) a, (0 + r, (t>) sj (t') , p (o) 

n »=1 
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where 



t ti 

dt'J dt"Ti E < 
o o 



k 



+ 



( rj (orj (O a, (f)aj(Op (o) + h.c. 



^+r, (or}(O5T(O5 J -(t ff )p(0) + h.c. 
' rj(or i (r)5 < (05t(r)p(o) + h.c. 
v +r, (or i (O5l(O5j(* ff )p(0) + h.c. 

m=i f r] (f ) a, (f ) p (0) rj (*") a, (*") + h.c. 

+r 4 (?) sj (0 p (o) rj (r) a, (t") + h.c. 



E 



rut') ^ (f) Ho) r J (t")^(t") + h.c. 



V+r i (05j(f)p(o)r J -(r)aJ(r)+h.c 



w\ ~t 



a? (t) — 6 1 aj, 



(8) 



and h.c. stands for Hermitian conjugate. 

If at t — the system is prepared in state ps (0) and the environment is 
admitted to be in thermal equilibrium, 



p(0)= Ps (0)®p E (0)=p(0) 



(9) 



where 



k n=0 



Pe (0) = - JJ exp (-(3hw k c\c k ) = II ( J 

oo 

Z = [[^exp (-(3f\w k n) , 



1 



it 

fik = 



+ n k ) \ 1 + n k 
1 

exp (/3huj k ) — 1 ' 



(10) 



/3 = (k B T)~ , k : B is Boltzmann constant, T the absolute temperature, Z the 
partition function and n k the mean number of excitations of frequency u k [12]. 
Taking the limit of zero temperature, it is easy to see that 



p s (t) - p s (0) = £) | dr [ f £ W* fe e^] e^V^-^*' (a,p 5 (0) aj - p s 

i,j=l o ^ fe ' 

+h.c, 

We used Fock states basis to take the traces over the environment and used 
for the integrals r = t' — t". 



Notice that J2k <^ikCi* k e lWkT decays very fast when r grows. Considering a huge 
number of environment modes, with frequencies spread over a large interval of 
values and a weak coupling to each mode (small a ik ), the phases of the e tuJkT 
make the sums practically null when r is not close to zero. If r c is the time 
when J2k a ikO:* k e' luJkT have appreciable values, we can, for t ^> r c , modify the 
integration limits without changing substantially the result. Thus we obtain 

p s (t) - p s (0) = £ { (fcy + iAij) (a,p 5 (0) at - p s (0) a}*) (j^dt'e^~^') } + h.c, (12) 

i,.7=l 



where the real numbers fc^- and Ajj are given by 

Tc 

% + zA^ = ]T a ifc «; fc y" dre^*-"^. (13) 

fc o 

Differentiating both sides of eq. (12) and iterating leads to 

d . _ * { (^ + iAy) (a,p 5 (i) a] - p s (t) a]a.) e^~<*)* 1 

= E S V , 7 +h.c..(14) 

dt id=i{ +0 2 (k t ,+tA t:i ) 

The terms 2 (kij + iA^-) are similar to the ones that appear if we consider 

O 4 (Hint) terms. Since we disregard O 3 (H in t) , we will not take them in to 
account. 

With the help of 

j t ~p s (f) = e*** |l [H 5 , p s (t)} + j t p s (t)| e"* 1 ** (15) 

we return to Schrodinger picture and write the master equation 

f t p s (t) = C Ps (t) , (16) 

where the Liouvillian superoperator (i.e.: operator which acts on operators) 
can be decomposed in three parts, 

C = C 1 + C 2 + C l2 , (17) 



and the parts are given by 
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£i 
£12 



a}ai,» 



a^a 2 , 



fcn (2ai • a{ — •a|a 1 — a]ai«) + i (A n — ui) 
k 22 (2a 2 • a| - »a^a 2 - aJ,a 2 «) + i (A 22 - uj 2 ) 
k 12 (a-i • a!, + a 2 • a| - »a|ai - a|a 2 «) + 
k 2 \ (a 2 • a{ + ai • a| - «ala 2 - a|ai«) + 

^ (ai • ajj, - a 2 • a| - «a|ai + a|a 2 «) + 
^ (a 2 • a| - ai • a| - «aja 2 + a£ai«) + 



A i2 


- A 2 i 




2 


A21 


- A 12 




2 


A12 


+ A 21 


2 



a x a 2 + a 2 a!, > 



(18) 



A word about notation is in order: we use the conventional notation for su- 
peroperators [13] where the dot (•) indicates the place to be occupied by the 
operator on which the superoperator is acting. 

Now we consider the environment large enough to take the limit of dense spec- 
trum. In this sense we define a density of modes D (uj) for the environment, 
where D (uj) duj represents the number of modes in the frequency interval be- 
tween u and 00 + duo. This allows the conversion of sums over modes (J2k) into 
integrals over frequencies (/ D(u)du), with the corresponding change of 
for a (uj) and /3k for (3 (uj), leading to 



kij + zAjj = J duD (uj) ctj (a;) a* (uj) £ 

Tc 

Q = Jdre^-^. 



(19) 



Using eqs. (3) we may write £ in terms of the operators related to the original 
laboratory modes A and B: 



a^a, 
Vb, 



£ — £a + £b + £inti 

£a = K a (2a • a f - «a f a - a f a«) + i (A aa - uj a ) 
£ B = h b (2b • b f - «b f b - b f b«) + i {A bb - u b ) 
£int = k ab (a • b 1 " + b • a 1 " — wb^a — a*b«) + 
k ba (b • a f + a • b f - «a f b - b f a«) + 

(a • b f - b • a 1 " - -b^ + a f b«) + 
(b • a f - a • b 1 " - «a f b + b f a«) + 



A 



ab 



*ba 



2 

A 6a - A 



ab 
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A ab + A ba 

2 

Here the coefficients are 



b+a + a+b,.!. (20) 



k aa + iA aa = j duD (u) {a (u) a* (u)r] + a (w) (3* (u) v} , (21) 

k ab + iA ab = J duD (to) {a (to) a* (u) v + a (to) (3* (u) fi} , (22) 

k ba + iA ba = j duD (u) {{3 (u) a* (lo) V + P M (3* (u) u} , (23) 

k bb + iA bb = J dujD (u) {(3 (u) a* (w)v + P (u) (3* (u) ^} , (24) 
rj = & cos 2 9 + ^2 sin 2 9, /i = £ 2 cos 2 9 + ^ sin 2 9, v = - ^ 2 ) sin 9 c^S) 

It is important to note that all superoperators appearing in eqs. (20) form a 
closed Lie algebra (see table 1). 

Ca and Cb are formally similar to the Liouvillians of independent modes 
evolving subjected to environment and Cm has the form of an "interaction 
Liouvillian" . However the constants k aa and A aa appearing in Ca depend 
on the coupling constants to the environment modes of both modes, a (to) 
and (3 (to) (and analogously for k bb and A bb ). The a (ca) (3* (u>) and (3 (co) a* (to) 
terms appearing in eqs. (21) and (24) are related to the return from the normal 
modes when we write C in terms of operators of the original laboratory modes. 
Notice that (eqs. (19)) ku and An [i — 1 or 2) depend only on jctj {uj)\ 2 - In fact, 
it is more natural to see the whole system through its normal modes, and not to 
look at each cavity mode. This is what the environment does: due to the form 
of £i and £2, the values of a (u), (3 (u) and D (cu) that effectively contribute to 
the coefficients of eqs. (20) are the ones with frequencies around the normal 
modes frequencies u>\ and u; 2 , and not around the original laboratory modes 
frequencies u a and u b . 

The physical interpretation for the coefficients A xy and k xy (x, y = a or 
b) is straightforward if x — y: A aa and A bb concern the system's unitary 
evolution, changing the oscillation frequencies of the cavities (the Lamb shift); 
k aa and k bb are dissipation constants. For x 7^ y the coefficients are related to 
a communication channel between the cavities mediated by the environment, 
playing unitary and non unitary roles in the dynamics (notice that k ab + iA ab 
and k ba +iA ba are not necessarily equal nor complex conjugate). As will become 
clear in sections (3) and (5), the values of k ab and k ba are very important for 
decoherence control. In eqs. (22) and (23) we find terms with a (cu) (3* (u>) and 
(3 (oj) a* (uj) which integrals tend to be null if the phases of a (u) and (3 (u) 
are not correlated. It is probably the most common case for fields in cavities. 
We must emphasize that analogous situation were studied in refs. [14] and 
[15] where the authors discuss the problem of two electric dipoles interacting 
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with the vacuum of electromagnetic field. The situation is formally similar, 
just changing bosonic operators for spin operators. However, as the dipoles 
are the only possibility of coupling, that system naturally exhibit this phase 
relation, and one can show [16,17] that this system is equivalent to the dipoles 
only coupled to one dissipative field mode. In the case of cavities, although it 
is desirable the phase relation, even in the uncorrelated case k a b and kb a may 
have appreciable values if v 7^ (eq. (25)), which may be achievable if, for 
example, the density of modes around lo\ have a different behavior than around 
uj 2 . It is another effect that appear when we write eq. (20), going back to the 
original laboratory modes. Since the environment here is composed mainly by 
phonons in the mirrors of the cavities and by electromagnetic modes in the 
space around the cavities, it is not impossible to manipulate the density of 
modes. 



3 General solution 



Before going over to specific examples, we present in detail the general solution 
of the master equation (20). 

The master equation (20) can be solved using Witschel's technique [18] and 
the parameter derivation technique [19]. This technique allows one to determine 
coefficients 5i such that the following identity is valid 

e (7iAi+7 2 A 2 H h-y n A n )t _ e <5i(t)Ai e <5 2 (t)A 2 . . . g i5„(<)A„ ^26) 



where the Aj's are superoperators forming a closed Lie algebra and t is a 
parameter (in our case, time). The parameter derivation technique consists of 
the following procedure: 

(1) Derive both sides of eq. (26) with respect to t and get 
(j2 TiAi) exp (jr lt A t t) = 5, (t) A, f[ e s 



\i=l / \i=l / i=l 



+5 2 (t) e Jl(t)Al A 2 I]e i,(t)Al 

i=2 

+ ••• 

n-1 

+5„ (t) J] e 5l(t)Al A n e 5n(t)A " . (27) 
i=i 

(2) Use the similarity transformation 

e xB Ae -xB = e xp,.] Aj (28) 
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and the linear independence of {Aj} in order to obtain differential equa- 
tions for the parameters Si (t). 
(3) Solve the c-numbers differential equations and obtain the factorized evo- 
lution superoperator written in the right hand side of eq. (26). 

As an example of step 2 above, let us take the second term in the r.h.s. of eq. 
(27) 



e «i(*)Ai A = e «i(t)Ai A -&(t)Ai e «i(t)Ai 
= e «i(t)[Ai,«] A2e «i(t)Ai 

= F 2 (5 1 (t),{A i })e s ^ A \ (29) 

Analogously one can define a similar operation for all the other terms and 
define 

F^F^S, (t),S 2 (t),{A i }),..-,F n = F 3 (5 1 (t),S 2 (*),•• -A-i (t) , {Aj}) , (30) 
and write 

(E exp (j2 -YiAitJ = (^i A i + 6 2 F 2 + ■ ■ ■ + 5 n F n ) exp (j2 7*A^ (31) 

or, equivalent ly, 

n 

£ 7iA 4 = <JiAi + 5 2 F 2 + • • • + 5 n F n . (32) 
i=i 

Comparing coefficients of each Aj (due to linear independence) one then ob- 
tains a system of coupled differential equations for the (Vs. 

In our case we have 
Ps (t)=e ct Ps (0) 

_ e h! (t)a.at e ft 2 (t)b«bt e z t (i)a.bt e «(t)b«at g n ( (t).atb e n(t)bta« 

e m 2 (t)btb» e p2(t)»b t b e mi(i)ata» e pi(t).ata e (j(t)atb» e g i (i)»bt a ^ ^ ^33^ 

Using the method just described we get 



^aa - k aa = m x (t) -n(t)q (t) e m ^~ m ^ , 
-iQ bb - k bb = m 2 (t) + n(t)q (t) e «*i (*)-««(*), 
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- t n ab -k ab = q(t) e "»i(*)-»*W 

-tQ ba -k ba = n (t) + n (t) (mi (*) - m 2 (t)) - n (t) 2 g (t) e ™i W 
+ifi aa - k aa = p, (t) - n x (t) q x (t) ePiW-wW, 
+iQ 66 - A: 66 = p 2 (t) + m {t) gi (t) ePiW-^W 
+ifi a6 -A; a6 = g,(t) ^W-wW, 

+ift 6a - A; fca = (t) + m it) (p! (*) - p 2 (*)) - (t) 2 g, (t) ewW-wW , 
2& aa = z (t) {-iVt ba - k ba ) 

+ Zl (t) (+in ba - k ba ) + h (t) (-2k aa ) + h (t) , 

2k bb = z (t) (+iQab ~ Kb) 

+ Zl (t) (-in ab - k ab ) + h 2 (t) (-2k bb ) + h 2 (t) , 

+i (Q ab - n ba ) + k ba + k ab = z (t) (+i£l aa - i^bb - k aa - k bb ) 

+hi (t) (-iftab - Kb) + h 2 (t) (+iQ ba - k ba ) +z(t), 
+i (Q ba - Qab) + K b + k ba = Zl (t) (-iVL aa + iVt bb - k aa - k bb ) 

+ K (t) (+in ab - Kb) + h 2 (t) {-iVtba - k ba ) + Zi (fjj34) 



where 



^aa — LO a — A aa , Qbb — Ub — Abb, fl a b — g — A ab , ^ba ~ 9 ~ A ba - (35) 



The solution reads 



n(t) 



h (t) 

h 2 (t) 
z{t) 



h (t) 



h it) 



h (t) ' * vv h it) 
. (t) , e m 

\f2(t)\ 2 +\l 2 (t)\ 2 



e m l(t ) = e -Rt fi (t) ^ e m 2 (t) = e - 2 Rt e - mi (t)^ 



\fl(t)\ 2 +\h (t)| 2 



1, 



1, 



-h (t) n (t) - 1* 2 w h (t) 

(3 ^r^m L 



Zl (t)=z* (t), 

n, (*) = (" (*))*, qi(t) = (q(t))*, 
p 2 (t) = (m 2 (t))\ Pl (t) = (mi (*))* , 



(36) 



where 



/i(*) 
h(t) 



1 

27 
1 

27 



r + c) e r + (r — c) e 



-rt 



(r - c) e + (r + c) e 



-rt 
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h (t) = — (ittab + Kb) (e 
h (t) = y (i^ba + ha) (e 



-rt 



-rt 




(37) 



and 



c 



^66 ^aa .(^66 ^oa) 



— +! 2 



r 



c 2 + {iVL ha + A; ba ) (zfU + k ab ), 



R 



kbb + &aa . ^ao + ^fcfe 



I Z =K m + 



(38) 



in 



The above equations look complicated. However there are some ingredients 
in them which become very illuminating from the physical point of view. Let 
us then discuss a few of them in a general context and apply the results in 
specific situations in the next two sections. 

It is probably sufficient to discuss the equation 



since it carries the essential ingredients of dynamical evolution. The point we 
want to stress is that by controlling some parameters, one can have distinct 
regimes of dissipation and decoherence in this model. The factor e" Rt in eq. 
(39) gives the expected decay behavior, with the average dissipation rate k m . 
The term e~ rt contains two contributions: a term which measures to which 
degree the two cavities are different, c, which adds to the other term coming 
from the dynamical coupling of the modes through the environment (k ao — iA ao 
and kba — i^ba) an d directly, i.e.: "through an antenna", (ig). One interesting 
effect is that for some parameter choices, this e rt term can cancel the deleteri- 
ous effects of the overall factor e~ Rt . For example, take Q aa = Qb and g — 0, 
i.e.: two degenerated modes. Note that if (k au — iA a b) (kj, a — iA^a) = k aa k b, 
we get r = k m , and, in this term, dissipation can be partially suppressed. 
Decoherence control is a little more intrincated matter, since it also depends 
on initial states as well. This will be discussed shortly. 




(39) 
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4 Two weakly coupled cavities with initial coherent states: manip- 
ulating dissipation 

In order to analyze specific cases, let us start with the simple situation in 
which coherent states are initially fed into the cavities: 

Ps(0) = \v a ) (v a \ <g> \v b ) (v b \ . (40) 

Applying the results of the previous section we get for the time evolution of 
this state 

p s (t) = \v a (t)) (v a (t) \®\v b (t)) (v b (t) | , (41) 
where 

Va (t) = e-( kM ^ {v a h (t) + V b h (*)} , 

v b (t) = e -(*™+^)* { V J 2 (t) + v b f 2 (*)} . (42) 

The first point to stress is that the known result of dissipative zero temperature 
time evolution of coherent states linearly coupled to the environment gives 
coherent states [20] applies here. Another point is that the auxiliary functions 
fi and li are oscillatory, generalizing the well known behavior of cosine and 
sine in simple cases like two identical isolated modes. Finally, as it should be, 
the symmetric character of the expressions: each mode fed the other with a 
strength depending on the coupling ifl ab + k ab and iVL ba + k ba . 

In order to gain further insight let us make some simplifying but realistic 
assumptions. Consider two modes in different cavities (A and B) with the same 
frequency, uo a = uj b = uj. The cavities are distant in the scale of the wavelength 
of the relevant mode, what leads to uncorrelated a (u) and (3 (u) (see section 
(6)) and then / D (u) a (u) (3* (u) & = (i — 1 and 2). If we have also D (cu), 
\cx(u)\ and slow varying around u\ and u 2 , the lamb shifts will be 

negligible: A aa = A bb = 0. Consider now very different dissipation rates, 
k bb ^> k aa (two geometrically identical cavities made of different materials, 
e.g.: one superconductor and the other metallic), and weak direct coupling 
in the sense g Ak (Ak = k bb — k aa ). Small g implies ui\ ~ u 2 — oo. Thus 
v ~ and, since we already have / D (u) a (to) j3* (to) ^ = 0, the coupling 
through the environment practically vanishes: k ab — iA ab = k ba — iA ba = 0. 
These conditions lead to 

c _ A A; _ 2g 2 
r ~ VAF - Ag ~ + AF' 
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(43) 



and consequently to 



Vb ( t ) = e -*-« {e- fe -^ b + 2| (e"*-' - e"M) Ua J , (44) 



where 



A:+-A; aa + — , 



^2 



k - = hb ~^k- (45) 

The natural decay rate for cavity A (B) is approximately k aa (k bb ). Under the 
assumption we made (k bb ^> k aa ), if they were uncoupled cavity B mode would 
decay much faster than cavity A. The effect of the coupling is to alter the decay 
rate of both cavities, leaving in cavity B a component (which can be made 
large by controlling v a ) which will decay with approximately the rate of cavity 
A. This result resembles the classical situation of pendulum synchronization, 
in which the center of mass motion is superdamped and the worst pendulum 
becomes phase locked to the best one. 



5 Superposition states: controlling coherence 

We next model a feasible experimental situation where one constructs entan- 
gled states between two cavities [4]. We consider an initially pure entangled 
state 

|V(0)> = W3(| W ,0)-|0,«;)), (46) 

where w is an index for coherent state, N — (2 — 2 exp |iu| 2 )) and clearly 
the density operator is 

p(o) = |v(o)> (0)1 - 

The time evolution here is given by 

p(t)=N{(\a 1 (t), e 2 (f)) - h (t) , a 2 (t))) ({a, (t) , e 2 (t) \ - e"* (e 1 (t) ,<J 2 (t)\)- 
'|<0H| 2 -|(ei (t)\a 1 (t))(e 2 (t)\a 2 (t))\e- i ^ 



(ei (t)k (t))(e 2 (t)\a 2 (t)) 



1^ (t))(e 1 (t)|<g>|e 2 (t))(a 2 (t)|+h.c.(4 
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where 



a, (t) 


= e 






a 2 (t) 


= e 


-k m t 




ei (*) 


= e 


-kmt g~ 




e 2 (t) 


= e 


~kmt g" 





are also indices for coherent states. Note that the dynamics of the non-diagonal 
elements is dictated by the difference between the "distances" between the 
components of the initial state, |(0|w)| 2 , and a similar quantity for the "time 
evolved" states in both cavities. The condition for coherence preservation is 
that 

\(0\w)\ 2 = (e 1 (t)\a 1 (t))(a 2 (t)\e 2 (t)), (49) 



i.e.: there is time evolution, but it is such that the scalar product of the 
evolved states is preserved. An extreme case would be unitary evolution, but 
we enforce that we only need this "unitary-like behavior" with respect to the 
vector states evolved in the initial state (46). 

A simple way to study state's purity is its linear entropy 

5{p)=Tr{p-p 2 }. (50) 



For the state of eq. (47) one gets 

5(t) _ (y 2 (t)-l)(x*-yZ(t)) 
2y*(t)(l-x) 2 ' 



where x = |(0|-u;)| 2 and y (t) = |(ei (t) \a\ (t))\ \(a 2 (t) \e 2 (t))\. This latter quan- 
tity generalizes the well known result that the time scale of coherence loss in 
a single cavity depends on the "distance" of the states which constitute the 
initial superposition. Time dependent "distances" such as |(ej (t) \<Ji {t))\ 2 ap- 
pear in the analytical expression for 5 (t). Note that S (t) — can be obtained 
by setting k aa = k bb = k ab = k ba , u a = u b and A aa = A bb = A ab = A ba = 0. 



6 Experimental (im) possibilities 



It is clear that any "unexpected" result in the problem of two modes inter- 
acting with one reservoir comes from the cross terms k ab , k ba , A ab and A ba . 
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In fact, as usual in quantum mechanics, one is dealing with an interference 
phenomena. In ordinary two ways interference phenomena (e.g.: double slit), 
in order to achieve large visibility, one tries to balance the interferometer and 
to guarantee coherence between the two ways. Here, the interference is due 
to dissipation rates. Balance means that each reservoir mode interacts with 
both oscillators with the same strength. In a two ways interferometer, coher- 
ence means a fixed phase relation between the two ways. The same situation 
applies to coupling constants here: complete "coherence" (we will use "coher- 
ence" whenever we use the word coherence with this meaning; there must be 
no confusion with the term decoherence) is achieved only when the phase of 
^ is independent of k. When the modes are directly coupled (g ^ 0), another 
effect can also originate cross terms, associated e.g. to variations of the density 
of modes between the two normal modes frequencies [21]. 

We begin by considering two field modes of distinct cavities. In principle, 
the reservoir contains modes of all frequencies. It is not bad to consider that 
the large wavelength modes interact with both cavities' modes in a "coherent" 
way (something similar to the dipole approximation [22]). However, the modes 
close to resonance are potentially the most important ones, as is emphasized 
by eq. (13). If on top of that the most important reservoir modes are also 
field modes, this means that the most important wavelengths are close to the 
wavelength of the cavities' modes. This suggests difficulties to implement this 
situation with distant cavities in the mode wavelength scale. Finally, the short 
wavelength modes of the reservoir, even without very large contributions to 
individual decay processes, generally will tend to destroy this "coherence" 
condition. As long as no directly coupling is imposed, one doesn't need to 
consider other effects. 

The above argumentation suggests the use of two modes of one and the same 
cavity. Now, if one thinks in terms of wavelengths, the most important reser- 
voir modes may possibly interact with the cavity modes in a "coherent" way. 
However, two (almost) degenerated modes form a two level system, and one 
can think about them as two polarizations of one spatial mode. The "coherent" 
interaction with the reservoir then defines a superdamped and a subdamped 
polarizations, in a way similar to super and subradiant modes of an n atoms 
maser [23], and in the limit situation, the subdamped mode would become a 
decoherence (and dissipation) free subspace (or mode). But what if we change 
our description from the previous polarizations to the super and subdamped 
ones? Then we would be describing a system where only one polarization is 
coupled to the reservoir, while the other is (almost) free. This is not the case, 
in general, but it is the case for some systems. One good example is the ex- 
periment of ref. [24], where DFS principle is addressed. A classical analog is 
given by two coupled harmonic oscillators, where center of mass mode usually 
couples strength with the environment and decays (much) faster then the rel- 
ative motion. Another favorable situation is the dipole situation [14,15] where 



16 



one has a first principles model for the coupling constants, and "coherence" 
can be naturally achieved manipulating the dipoles. 

Obviously an explicit mechanism regarding the geometry of the cavities is still 
missing. It will be introduced and compared to experiment in a forthcoming 
publication [25]. 

In short, we studied the effects of correlations among coupling constants of 
each of two modes with a common reservoir, showing that they can smoothly 
vary from two independent damped systems, to one superdamped and one free 
modes, in case of perfect correlations and/or ubiquitous environment density 
of modes. Examples of state evolution were given and some experimental con- 
siderations were made, and we believe that some light was shed on decoherence 
and dissipation processes for composed systems in cavity QED when the pro- 
posed decoherence mechanism seems to work well [11]. 
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Table 1 



Commutation relations [xi,Xj] f° r the bosonic superoperators (xi and Xj are given 
by the line and the column respectively). 
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